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ABSTRACT 

Based on a deep imaging survey, we present the first homogeneous star formation history (SFH) of 
the Fornax dwarf spheroidal (dSph) galaxy. We have obtained two-filter photometry to a depth of 
-B ~ 23 over the entire surface of Fornax, the brightest dSph associated with the Milky Way, and derived 
its SFH using a CMD-fitting technique. We show that Fornax has produced the most complex star 
formation and chemical enrichment histories of all the Milky Way dSphs. This system has supported 
multiple epochs of star formation. A significant number of stars were formed in the early Universe, 
however the most dominant population are the intermediate age stars. This includes a strong burst of 
star formation approximately 3 — 4 Gyr ago. Significant population gradients are also evident. Similar 
to other dSphs, we have found that recent star formation was concentrated towards the centre of the 
system. Furthermore, we show that the central region harboured a faster rate of chemical enrichment 
than the outer parts of Fornax. At the centre, the ancient stars (age > 10 Gyr) display a mean metallicity 
of [Fe/H] ~ —1-4, with evidence for three peaks in the metallicity distribution. Overall, enrichment in 
Fornax has been highly efficient: the most recent star formation burst has produced stars with close to 
solar metallicity. Our results support a scenario in which Fornax experienced an early phase of rapid 
chemical enrichment producing a wide range of abundances. Star formation gradually decreased until 
~4 Gyr ago, when Fornax experienced a sudden burst of strong star formation activity accompanied 
by substantial chemical enrichment. Weaker star forming events followed, and we have found tentative 
evidence for stars with ages less than 100 Myr. 

Subject headings: galaxies: individual (Fornax dSph) — stars: formation — Local Group 



1. INTRODUCTION 

The dwarf spheroidal galaxies (dSphs) are the least lu- 
minous galaxies known. They display remarkably high 
mass-to-light ratios (~100 to 1000), and the stars in each 
system are known to reside at the centre of a massive dark 
halo (Mvir 10^ - 10^ M(7)) which ex tends far beyond the 
obser ved limiting radii (jWalker et al.|[2007. : Simon fc Gehal 
l2007t ). In terms of stellar population, these systems can 
have surprisingly complex star formation histories (SFHs). 
All d Sphs contain a population of ancient stars ( Held et al.l 
120001 ) ■ however some (such as Fornax) have been able to 
maintain multiple epochs of star formation and chemical 
enrichment over a Hubble time. These systems are rel- 
atively simple environments compared to larger galaxies, 
and are therefore a starting point in the study of star for- 
mation and enrichment. Simulations suggest a cyclical 
process, in which the gas collapses to form stars, is then 
chemically enriched and blown out by pockets of massive 
star formation, and then collapses again to repeat the cy- 
cle. iSalvadori et al.. (,2008. ) propose a time frame of ~250 
Myr for a single cycle. 

There are, however, open questions regarding star for- 
mation in dwarf galaxies. Population gradients sug- 
gest that the most recent bursts of star formation in 
each dSph were conc entrated towards the object's centre 
(|Harbeck et al.ll2001h . Also, the Milky Way Halo contains 
a population of extremely metal-poor stars, with abun- 
danc es of [Fe/H] < -5 (Christhcb et al. 2002; Frc bel et al.l 
l2005f ). whereas the dSphs do not co ntain a populatio n 

with metalhcities below [Fe/H] 3 (jHelmi et al.ll2006 D. 

This suggests that the gas sourcing the first generation 



of stars in these systems was pre-enriched. Moreover, 
iGrebel fc Gallagher! (|2004D have argued that the variety of 
SFHs in dSphs is not the result of reionization, and hence 
'local processes' have influenced star formation in each ob- 
ject. These includ e the regulation of gas dynamics due to 
internal feedback (jPekel fc Silkl [19861 ) while ram pressure 
stripping and tidal interaction with th e Milky Way are also 
important factors (Ma ver et al.ll2006l ). 

With an in tegrated absolute magnitude of My = —13.1 
(|Mateolll998l) . Fornax is the brightest dSph associated with 
the Milky Way (excluding the tidally dirsupting Sagit- 
tarius system). It lies at a distance of 138 kpc, and 
proper motion measurements indicate its orbit around 
the Milky Way is roughly circular with an eccentricity 
of e = 0.13j |nnT and an orbital period of 3.21q7 Gyr 
(jPiatek et al.ll2007t ). IWalker et all (|2006[ ) have completed 
a large kinematic survey of this system, collecting radial 
velocities for 206 member stars over the entire surface of 
Fornax. They measured a flat velocity dispersion profile, 
indicating a significant dark component, and find a mass- 
to-light ratio oi M/Ly ^ 15 within 1.5 kpc (approximately 
half the tidal radius). 

Compared to other dSphs, the star formation his- 
tory of Fornax is unusually comple x. All dSphs con- 
tain some number of ancient stars (jCrebel fc Gallaghed 
[2004") , however early studies of Fornax revealed an ex- 
tended giant branch, including carbo n stars, suggesting: 
a strong intermediate-age componen t (jPemers fc Kunkel 
[T979I: lAaronson fc Mouldl[T980L [19851) . Additionally, For- 
nax contains a sign ificant young stellar com ponent. Origi- 
nally discovered bv lBeauchamp et al.l (|l995f ). the photom- 
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etry of lStetson et~all (|l998l ) and lSaviane et all (|2000f ) later 
identified main sequence stars with age s as yo ung as 100- 
200 Myr. The analysis of iGallart et al.l (|2005[ ) indicates a 
burst of star formation in the centre of Fornax 1-2 Gyr 
ago, which has continued almost to the present day. In- 
deed, more than half the s tars on the RGB are thought 
to be younger than 4 Gyr (jBattagha et al.ll2006D . Of all 
the dSphs, Fornax has experienced the most recent star 
formation. 

Fornax also contains population gradients, in which the 
young stars ar e preferentially located towards the centre 
(jStetson et al.| [l998: Battaglia et al. 2006). T his property 
is common in the dSph population ( Harbeck et al.|[200l[ ) 
and suggests that the gas required for subsequent star for- 
mation episodes was more successfully retained in the core 
of the dark halo than in the outer regions. In Fornax, the 
young component is not aligned with the main body and is 
highly structured, including a shell-like fea ture which may 
indicate an accretion event ~2 Gyr ago (jColeman et al.l 
12004 lOlszewski et al.l[2006h . 

In addition to an age spread, the star s in this sys- 
tem c o ver a signifi c ant ra nge in metallicity. Tolstoy et al.l 
(|200H ). lPont et al.l (120041 ) and lBattagha et all (|2006[ ) have 
examined the chemical abundances of Fornax red giants 
using spectra of the Ca II triplet features, and found a 
metallicity range of —2.5 < [Fe/H] < 0.0. The data set of 
Battaglia et al. contained enough Fornax members (562) 
such that stellar metallicities could be accurately related 
to kinematics. They found that the metal-rich stars have 
a colder velocity dispersion, and the metal-poor compo- 
nent shows signs of non-equilibrium kinematics towards 
the centre of Fornax (r < 2rc). Additionally, high reso- 
lution spectra of 81 red giants in the centre of Fornax in- 
dicate that s-process elements are unusually strong, hence 
stellar winds (such as those from AGB stars) have domi- 
nated th£_cheinical enrichment of Fornax in the last 2-4 
Gyr (lLetartel[2007[) . 

Despite the progress described above, Fornax has been 
lacking a homogeneous study of its SFH over the en- 
tire surface of t he sy stem. The HST photometry of 
iBuonanno et al.l (|1999D confirmed the presence of young 
stars at the centre of Fornax and included evidence of 
seperate bursts of star formation. B y combining these re - 
sults with Ca II triplet metallicities, Tolstoy ct al. (20011) 
created schematic star formation and chemic al enrichment 
histor ies. The preliminary SFH produced bv lGallart et al.l 
(|2005D was based on VLT/FORSl photometry (depth 
/ ~ 24.5) also located at the centre of Fornax, and an- 
other field approximately one core radius from the centre. 
Young stars were present in both fields, however the au- 
thors noted a significant difference in SFH between the 
two. In summary, although the general trend of star for- 
mation at the centre of Fornax is known, the aggregate 
history is yet to be determined. 

Hence, we present the first results of a deep, homoge- 
neous photometric data set over the entire surface of For- 
nax. We have extracted the SFH from this photometry 
using the CMD-fitting techniques developed by Dol phinl 
(|2002D . Not only has this allowed us to derive a global 
SFH for this system, we have also examined the SFH as 
a function of position to search for population gradients. 
Our data have a limiting magnitude of i? ~ 23.0, thus we 



are sensitive to main sequence stars in Fornax with an age 
of 3 Gyr and less. Also, the red giant branch, red clump 
and horizontal branch stars allowed us to track the ages 
and metallicities of the populations with ages > 3 Gyr. 
This is the first complete SFH of this system derived from 
deep photometry. 



2. THE SURVEY 

2.1. Data Reduction 

Images encompassing the surface of the Fornax dSph 
were obtained using the ESO/MPG 2.2m Telescope 
equipped with the Wide-Field Imager. This instrument 
provides a 34' x 33' field of view using a 4 x 2 mosaic of 
2048 X 4096 pix^ CCDs and a pixel resolution of 0.24" 
pix~^. The survey aim is to obtain photometry to 23rd 
magnitude over the entire body of Fornax. The first stage 
presented here contains Fornax itself and the outer shell 
noted by IColeman et al.] (|2005t) . Thus far, 21 pointings 
have been obtained, covering a sky area of 5.25 deg'^. An 
overlap region of ~4' between each field was chosen to 
ensure the photometric zeropoint was constant across the 
survey. A schematic diagram of the fields is shown in Fig. 
[TJ We obtained three 600s dithered exposures in both 
B and R for all fields, allowing us to reach magnitudes 
of B = 23.0 and R = 23.5 (50% completeness limits) 
in all fields. The images were taken during 11 nights in 
Oct /Nov 2006 in median seeing conditions of 1.4" (range 
0.9" — 2.1"). The data were reduced using standard proce- 
dures in the mscred package in IRAF: the overscan region 
and the bias frames were used to subtract the pedastal cur- 
rent from each science image, which were then trimmed. 
Twilight flat fields were combined to produce a master flat 
frame in B and R for each night, which were then used to 
flat-field the science images. 
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Fig. 1. — Schematic diagram detailing the layout of the fields 
curre ntly observed. The red ellipses represent the core and tidal 
radii UMateol IToM) . The dashed lines are the limit of our sur- 
vey and the dotted lines represent the overlap regions between 
fields. The shaded circles represent the five globular clusters as- 
sociated with Fornax and are labelled accordingly. These circles 
have radii of 1.5', approximately equal to the tidal radii of each 
cluster (Mackcy & Gilmorc 2003). The two arcs represent the shell- 
like features noted by Coleman ct al. (2004, 2005), and the dashed 
ellipses represent the boundaries of our radial bins. 

Three dithered images were taken in each field, and we 
used routines in the mscred package to co mbine the im - 
ages and remove gaps between the CCDs (IValdesI l2002l ) . 
An astrometric solution was constructed for all images by 
matching them to the first USNO CCD Astrograph Cat- 
alogue (UCACl), which has an average precision of 31 
mas in the magnitude range 8 < i? < 16 (jZacharias et al.l 
l200Clf ). The rms of our solution was <0.2" in all fields. 
The individual CCD images were then combined using the 
mscimage routine to produce three single (i.e. non-mosaic) 
images for each field in both filters. These three images 
were then median combined with the mscstack routine, 
which matches images based on their astrometry and re- 
moves the CCD gaps. Finally, the combined image for each 
field was then corrected for zeropoint gradients across the 
field of view using the mscskysub routine. 

2.2. Photometry 

Ph otometry was derived using DAOPHOT (iStetsonI 
Il987f ). We measured the background level in each field to 
calculate a standard deviation of the sky, a. Each image 
was then searched for all sources Aa above the background 
level, and aperture photometry was used to estimate their 
brightnesses. In a crowded field such as the centre of For- 
nax, the PSF-fitting technique in DAOPHOT provides a 
more accurate measure of stellar magnitudes compared to 
aperture photometry, as it allows the signal of adjacent 
stars to be disentangled. Hence, we examined the bright- 
est 60 stars in each image and used those with no apparent 
neighbours and a well-defined Gaussian shape to construct 
a master PSF for each science frame. This was then fitted 
to every source in the image using a fitting radius of 8 
pixels (1.9") or 1 — 2 half-width, half-maxima of the PSF 
(depending on seeing). 

To measure the completeness and photometric accuracy 
of our photometry, we performed artificial star tests on 
all science frames in both filters (i.e. 42 science images). 
We placed 1600 artificial stars in the image and attempted 
to recover them with DAOPHOT, where the photometric 
uncertainty was then determined as the dispersion of the 
returned magnitudes about the mean (that is, not the in- 
put). This was repeated for artificial stars at every 0.25 
magnitudes in the B and R frames. 

To ensure a constant zeropoint across the survey, we 
matched stars in the overlap regions between fields using 
their astrometry and measured the mean inter-field differ- 
ence in B and R. This i s the same technique used in our 
previous Fornax survey man et al.ll2005D . The inter- 
field corrections were accurate to ~0.02 mag in both B 
and R. As an example, the final match between the fields 
F6 and Fll is shown in Fig. [2l A final zeropoint cor- 
rection was made by matching our data to the catalogue 

1 http://www.ls.eso.org/lasilla/sciops/2p2/E2p2M/WFI/filters/ 



of lStetson et al.l (jlQQSD . which contains photometry of the 
core of Fornax in B and i? to a similar limiting magnitude 
as our survey. The R filter attached to the 2.2m/WFI is a 
standard Cousins filter, and the i?-band photometry was 
well matched between the two datasets. In contrast, the 
B filter (BB#B/123J]S0878) covers a larger wavelength 
range than the standard Johnson filteiQ, and therefore re- 
quires a colour correction to match standard B magni- 
tudes. The ESO website provides a correction for the B- 
band photometry using the (B — V) colour; we determined 
a {B — R) colour correction by comparing our data to the 
Stetson et al. catalogue, yielding the following result, 

B ^b + 0.20{b~ R), 
with a bootstrap error of 0.02 mag. The overall photomet- 
ric zeropoints are accurate to 0.03 mag. 
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Fig. 2. — A comparison of the photometry measured in fields 6 
and 11. We have removed all sources with a non-stellar sharpness 
and large photometric uncertainty (cja.^) > 0.1 mag) 

2.3. Colour Magnitude Diagram 

The colour magnitude diagram (CMD) of Fornax rein- 
forces the complex star formation history (SFH) of this 
object: it contains multiple stellar populations with a vast 
range in age and chemical abundance. A full description 
of the stellar population s at the centre of Fornax is given 
bv IStetson et all (| 199^ and ISaviane et all (|2000[ ). How- 
ever, the outer regions of Fornax are not well known, and 
in Fig. [3] we present the first deep CMDs for the entirety 
of Fornax. It has been theorised that For nax may con- 
tain s trong radial population gradients (e.g. lSaviane et all 
1 2j)0_C|) , hence we have divided the data set into the four 
elliptical regions shown in Fig.[TJ The regional boundaries 
are at radii of rc, 2rc, 3rc and rt, where rc = 13 .8' and 
Tt = 76.0' are the core and tidal radii listed bv iMateol 
()1998f ). The distribution of young stars in Fornax is not 
aligned with the system's major axis and it is also know n 
to contain strong asymmetries (e.g. IStetson et al.l Il99"8l ). 
however the majority of these stars are contained in the 
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core region (this will be shown in the next paper in this 
series; Coleman, in prep.) and hence the assumption of 
elliptical regions is not vital for this sub-population. 
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Fig. 3. — Colour-magnitude diagrams of the four Fornax regions. 
The upper four panels show the photometry for sources in Region 
1 (r- < rc), Region 2 (r'c < r < 2rc), Region 3 (2rc < r < 3rc) 
and Region 4 (3rc < r < rt). The photometry is 95% complete to 
a magnitude of B = 23.0, however this number improves slightly 
(~0.25 mag) in the outer regions due to less stellar crowding. We 
have removed all sources with a non-stellar sharpness and a large 
photometric uncertainty {(^(^b—R) > 0-3 mag). The four lower pan- 
els represent the 'signal' of each subset above the background. 

Region 1 (the core region) contains approximately 
28,000 sources. Moving outwards, another 41,000 stars 
were selected in Region 2, 31,000 in Region 3, and 51,000 
in Region 4. For convenience, we have also performed a 
background subtraction on each CMD, removing the field 
po pulation using the 'sig nal-to-noise' technique provided 
bv lGrillmair et all (|l995f ). A full description of this tech- 
nique as applied to the current data set will be given in 
the next paper in this series. In summary, we have di- 
vided the CMD into a grid of cells, counted the number of 
Fornax and field stars in each cell, and then used this to 
remove the field contamination (mostly foreground Milky 
Way stars). The results are shown in the lower panels of 
Fig.H 

The central CMD (upper left panel) shows all sources in 
Region 1, drawn from within the inner red ellipse shown in 
Fig- El Immediately visible is the red giant branch extend- 
ing downwards from B ~ 20, which contains old and inter- 
mediate age stars (age > 1 Gyr) and is thickened due to 
the age and metallicity spread. The most densely packed 
feature is the red clump (RC), centred at {B ~ R) — 1.3, 
B = 22. This is made up of core helium burning stars, 
and is essentially a young-to-intermediate age, metal-rich 
horizontal branch (HB). A hint of a HB extending towards 
the blue is also visible. We also note that the HB appears 
to extend redward from the RC, however these objects are 
artefacts of the observational dithering pattern: a CMD- 
selection indicated that they lie pre-dominantly in regions 



in which only a single frame of B-band data was available. 
This effect is most prominent in the HB star-rich central 
four fields, where a few of the HB stars display large pho- 
tometry errors in the B filter (and hence an artificial colour 
spread in the HB itself is created). Although it is a small 
artefact (it exists in less than 1% of the surveyed area), 
its effect on our star formation history results is discussed 
further in gSXSl 

Returning to our discussion of the central region CMD, 
a slight overdensity lies approximately 0.6 mag above the 
RC, identified as the c olour-magitude clumping at the 
start of the AGB cycle (jSaviane et al.l [20001. Below the 
RC lies the sub-giant branch, which is also thickened by 
the age and metallicity spread. Finally, the blue column of 
stars extending downwards from (B — R) = —0.2, B = 20 
is the young main sequence. iBeaucha mp et al. (1995*) dis- 
covered this feature in Fornax, and Savianc et al. (2000|) 
subsequently identified main sequence stars as young as 
200 Myr. Indeed, of all the Milky Way dSphs, Fornax has 
the most recent star formation. 

At first glance, the thick red giant and sub-giant 
branches are common to all four regions, indicating that 
the wide range of stellar ages and metallicities continues 
well beyond the core radius. However, the differences en- 
countered in the stellar population when moving outwards 
from the centre of Fornax are remarkable. We see a de- 
crease in the prominence of the young main sequence, pos- 
sibly indicating that recent star formation (i.e. less than 
4 Gyr ago) was preferentially located towards the centre 
of Fornax. Another clear difference between the four re- 
gions lies in the morphology of the HB. The red clump is 
present in all four regions, however the HB extends further 
into the blue region as we move outwards; it terminates at 
(B - i?) ~ -0.2 in the Region 4 CMD. This indicates that 
the outer regions of Fornax contain a significant (if not 
dominant) population of old, metal-poor stars. Further 
emphasis of this point is provided by the red giant branch 
(RGB): the lower panels of Fig. [3] show that as we move 
outwards, the RGB shifts towards the blue, thus indicat- 
ing a decrease in mean metallicity with increasing radius. 
Overall, a comparison of the four CMDs shown in Fig. 
[3] would suggest that later bursts of star formation and 
chemical enrichment were preferentially located towards 
the centre of Fornax. 

3. STAR FORMATION HISTORY 



Numerical fitting of CMDs allows a study of the 



SFH of a dwarf 


galaxy (e.g. ' 


Gallart et al. 


'1996; 


Tolstov & Saha 


1996; 


Aoaricio et al., 


1997; Dolphin 1997; 


Holtzman et al. 


199C 


Olsen 119991: Hernandez et al.ll2000l: 


Harris & Zaritskvl 200lh. Moreover, the large mosaic ore- 



sented in this paper provides an opportunity to study the 
spatial variation of the SFH within Fornax. To obtain a 
detailed picture of the SFH w e use the CMD-fitting soft- 
ware MATCH polphin|[200l . which applies maximum- 
likelihood methods to fit photometric data with simple 
model CMDs. By converting data and models to so-called 
Hess-diagrams (2-D histograms of t he stellar d ensity as 
a function of colour and magnitude; iHessI [192^ a direct 
pixel-by-pixel comparison is possible. The model CMDs 
are based on theoretical isochrones from iGirardi et al.l 
(|2002i) and include realistic photometric errors and com- 
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pleteness, which are obtained from the artificial star tests 
described earher. By determining the best-fitting linear 
combination of model CMDs for different age and metal- 
licity bins, the SFH and metallicity evolution are then con- 
strained. The accuracy of the recovered metallicities de- 
pends not only on the data quality, but also on the quality 
of the isochrones and the stella r evolu tion tracks on which 
they are based. Ide Jong et all ([2008) tested MATCH on 
a set of six globular clusters with varying metallicities us- 
ing i sochrones based on t he same stellar evolution tracks 
from lGirardi et all (|2002[ ). They show that the recovered 
metallicities are always within 0.2 dex of the spectroscopic 
values. In all results presented in this paper we therefore 
include an additional contribution to the metallicity un- 
certainties of 0.2 dex. 

Since there are variations in seeing and sky brightness 
between the different fields in the mosaic, the SFH fits are 
done separately for each field. Furthermore, the four ra- 
dial bins are treated separately to enable an analysis of 
the radial variation of the stellar populations in Fornax. 



3.1. CMD fitting method 

The main free parameters in CMD fitting are distance, 
age, metallicity and extinction, although the binary frac- 
tion and the assumed initial mass function (IMF) also play 
a role. To limit the number of free parameters to the age 
and metallicity, we assume reasonable priors on the others 
parameters. For all our fits we assume a binary fraction 
of 0.5 and a Salpeter IMF (|Salpeteil [19551 ). which to the 
CMD depth probed he re is practical l y equ al to, for ex- 
ample, a Kroupa IMF (|Kroupa et al.|[l993l ). During the 
past decade, several different studies have a ll found con- 
sistent distances to Fornax of 138 ± 4 kpc (jMateol[l998t 
lBersieij[2000t iRizzi et al.ll2007D , justifying a prior on the 
distance in our SFH fits. To account for the small un- 
certainty in the distance, we perform all fits for three 
fixed distances, namely 135, 138 a nd 141 kpc. Accord - 
ing to the dust extinction maps from lSchlegel et all (I1998D 
the foreground reddening towards Fornax varies between 
E{B — V) = 0.015 and 0.03 mag. Fits are performed for 
different combinations of foreground and internal extinc- 
tions, with the former fixed at E{B - V) = 0.015, 0.02, 
0.025 and 0.03 mag, and the latter at E{B-V) = 0.0, 0.1, 
0.2, 0.3 and 0.4 mag. 

For each individual fit, the distance and two kinds of 
extinction are fixed, while the star formation rates (SFRs) 
for the age and metallicity bins serve as the free parame- 
ters. Since the isochrones are spaced more evenly in log(t) 
than in t, the age bins are defined in log(t). They have bin- 
widths of A log t = 0.15, with the oldest bin corresponding 
to ages between 11 and 16 Gyr and the youngest to 10 
and 16 Myr, for a total of 21 age bins. The full metallicity 
range spanned by the isochrones, [Fe/H] = —2.4 to 0.0 is 
covered with 16 bins of 0.15 dex width. To account for 
contamination by foreground stars and faint background 
galaxies, a control field CMD is created from all star-like 
sources outside the limiting radius of Fornax, and used 
as an additional model CMD in the fits. All stars with 
-1 < B - i? < 3 and B < 23 and i? < 23 are fit, using 
Hess-diagram bin sizes of 0.16 in magnitude and 0.08 in 
colour. 
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Fig. 4. — CMD fit for the central radial bin part of survey field 
15. (a) Observed Hess diagram, with the grayscale indicating the 
number of stars per pixel, (b ) Hess diagram of the best-fit model, in- 
cluding theoretical models and a control field population, (c) Resid- 
uals after subtracting the model from the data, where darker pixels 
correspond to underprediction and lighter to overprediction of stars 
by the model, (d) Residuals scaled by Poisson a. 

Figs. |4] and [5] show the Hess diagrams of Regions 1 and 
3 recovered from Field 15 with the corresponding best- 
fit models and residuals. In general the fits are good, 
but some systematic problems arise in reproducing the ex- 
act shape of the HB/RC and its extent towards the red. 
This imperfect modeling of the HB /RC region is a general 
problem for theoretical isochrones, due to the complicated 
processes taking place during this phase in stellar evolu- 
tion. Newer isochrone sets than the one currrently used 
by MATCH (Girardi et al. 2002) should improve this situ- 
ation, but our current results are not strongly affected by 
these problems. However, as the details of the theoretical 
isochrones used do infiuence the exact values of metallic- 
ities and ages we recover, some extra uncertainty should 
be taken into account when interpreting the fit results. 
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Fig. 5. — Same as Fig.U but for the third radial bin. 

After running the 60 fits (3 distances, 4 foreground and 
5 internal extinction values) for each region, the values of 
the maximum-likelihood mea sure of goodness-of-fit Q (see 
lDolDhinll2002l : Ide Jong et al.|[2008.) of all fits are compared. 
The expected random (i.e. not due to actual SFH differ- 
ences) variance in this Q parameter, cr, when fitting a spe- 
cific CMD, is calculated using Poisson statistics and from 
Monte Carlo simulations using random drawings from the 
best model. All fits that have a value of Q within la of the 
best fit are considered as 'good' fits and used to construct 
the SFHs presented in the remainder of this paper. 

3.2. Results 

3.2.1. A Global View 

As a first pass, we examine the aggregate SFH of the 
Fornax system. Combining the SFHs obtained from all 
regions and fields gives the total SFH of Fornax, pre- 
sented graphically in Fig. [S]and in tabular form in Table[TJ 
Clearly, Fornax has a complex history and has been form- 
ing stars continuously over most of the age of the universe, 
with non-zero SFRs being found in most age bins. After 
the first stars formed more than 10 Gyr ago, a slow decline 
is detected until approximately 4 Gyr ago, when the star 
formation rate sharply increased for a period of ~1 Gyr. 
After this sudden intensified star formation episode, a very 
low level of star formation has continued until the present 
day. For comparison, we also show in Fig. [6] the schematic 
SFH derived bv Tolstov e t al. (2001) based on photometry 
and Ca ii triplet results. We have shifted the peak of their 
relative star formation rate to match our peak SFR. The 
results are well matched for the last few Gyr, however our 
CMD-fltting code is better able to extract detailed star 
formation rates from the early history of Fornax. We also 
show the total stellar mass formed during each bin in Table 
[Hand Fig.[7l From this, we calculate the total stellar mass 
formed in Fornax to be e.ltjjj X 10^ Mq, where we have 
integrated the Salpeter IMF down to a mass of O.15Af0. 
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Fig. 6. — Total star formation history {upper panel) and age- 
metallicity relation {lower panel) for Fornax. In both panels, the 
horizontal error bars indicate the widths of the age bins. The er- 
ror bars on the star formation rates are not independent; a higher 
SFR in one bin would be compensated by a higher SFR in adjacent 
bins. Metallicities are only shown for age bins with a SFR greater 
than 0.0005 Mq yr~^ and are SFR-weighted means. Both panels 
contain a comparison with previous work, represented by the dot- 
ted l ines. The uppe r panel contains the schematic SFH derived by 
Tolst ov et al.l 1 120011) scaled to match our peak star formation rate 
detection, and the lower panel contai ns our estimate of the e volution 
of the average metallicity derived bv IBattaglia et aTl II2006I ). 
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Fig. 7. — Stellar mass formed during each age bin, taken directly 
from the total star formation history (Fig. |BJ. 

Initially, the average metallicity of the stars shows little 
evolution. A large fraction of the ancient stars seems to 
have formed from pre-enriched g the mean metallic- 
ity in the oldest age bin (11 — 16 Gyr) is already [Fe/H] 
= —1.4. The metallicity only starts to increase signifi- 
cantly from this value at an age of around 5 Gyr. Overall, 
these results are in good agreement with previous studies 
of the SFH and metallicity of Fornax. We show in the 
lower panel of Fi g. [6l th e approximate age-metallicity re- 
lation derived bv IBattaglia et al.l (j2006f ) based on spectra 
of 562 RGB stars extending from the centre of Fornax to 
the tidal radius. This line represents our estimate of the 
evolution of the average iron abundance from Fig. 23 of 
[Battaglia et al. (2006), however we note that each age bin 
encompasses a wide range of stellar metallicities and thus 
this line should be regarded as a guide only. Although they 
are well within the uncertainties, our results are slightly 
more metal-rich than the spectroscopic results. Such a 
small offset in metallicity might be caused by slight inac- 
curacies in the isochrones, the photometric calibration, or 
the distance to Fornax. 

3.2.2. A Refined View 

Population gradie nts are a common feature of dSphs 
(jHarbeck et al.|[200l[ ). Generally, more recent star forma- 
tion events are more centrally concentreated and produce 
stars with a greater metallicity compared to the older stel- 
lar population. In Fig. [5] we show the total SFHs (left) and 
metallicity evolutions (right) for all four radial regions (e.g. 
for Region 1 the four SFHs of the portions in fields 10, 
11, 15, and 16 have been combined). Overplotted in each 
panel are the SFHs and metallicities for four individual re- 
gions, to give an idea of the field-to-field variation. Since 
the panels on the right only show the average metallicities 
and provide no insight into the metallicity distribution, 



Fornax dSph I: SFH 



7 



the combined fit results for the four radial regions with 
the complete age-metallicity grid are shown in Fig. [9l 
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Fig. 8. — Relative star formation histories {SFH,ie/t) and age- 
metallicity relations (AMR, right) per radial bin. From top to bot- 
tom the panels correspond to the central through to the outermost 
radial bin. As in Fig. \E\ horizontal error bars show the age bin 
widths and the errors on the SFRs of adjacent bins are not indepen- 
dent. The SFHs shown here are scaled to the average SFR in each 
radial bin. Overplotted on each panel, using colored lines, are the 
SFHs and AMRs of each radial bin in four individual survey fields; 
which fields are plotted is indicated in the panels on the right. 

There is a striking difference in the SFH when moving 
from Region 1 (top panels) to Region 4 (bottom panels). 
In the center of Fornax, the burst of star formation that 
occurred 3 — 4 Gyr ago stands out strongly. Moving out- 
wards, the proportion of stars produced during this burst 
decreases and in the outskirts are not found at all. Stars 
younger than 1 Gyr are also mostly found in the center, 
with hardly any significant star formation at these ages 
in the outermost radial bin. This kind of radial depen- 
dence, with young stars more centrally concentrated than 
old stars is a common characteristic of dS phs, and evi- 
dence for this was found b efore in Fornax (jStetson et al.l 
Il998t iBattaglia et al.ll2006t) . At all radii, the metallicity 
evolution follows the same pattern. It starts off at [Fe/H] 
~ — 1.5 and remains constant until ~5 Gyr ago. Then the 
metallicity starts to increase rapidly to —0.5, after which 
there is little indication of further enrichment. This rapid 
increase in metallicity coincides with the strong burst of 
star formation. 
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Fig. 9. — CMD fitting results for the four radial bins. From left 
to right the panels are for the central bin (labeled with a "1") to 
the outermost bin (labeled with a "4"). The grayscale indicates the 
star-formation rate in each age-metallicity bin, with a darker value 
corresponding to a higher value. 



Comparing the SFHs of individual fields within a radial 
bin (i.e. the coloured lines in Fig. [8]) shows no significant 
difference in the central bin. Since in the central part 
the dynamical timescales are shortest, the populations are 
well-mixed. In the outermost region there is also no sign 
of spatial variation, which is probably because there are 
mostly ancient stars which have had sufficient time to dif- 
fuse throughout the galaxy. In Regions 2 and 3, however, 
significant field-to-field variation is visible in the strength 
of the 4 Gyr burst. The detailed structural properties of 
the stellar populations in Fornax is the subject of the next 
paper in this series (Coleman, in prep.). 



3.2.3. Abundance Variations in the Ancient Stars 

Although Fornax is found to be relatively metal-rich, 
Fig. [5] shows that the spread in metallicities is significant. 
At all ages the metallicity spread is at least '^1 dex, and 
especially the oldest stars (>10 Gyr) show a very large 
spread, with metallicities ranging between [Fe/H] — —0.5 
to our metallicity cut-off at —2.4. Note that in all regions 
there is a non-zero SFR in the most metal-rich two age 
bins. This is caused by MATCH attempting to fit the red- 
ward extension of the HB (see Fig. [3]) with an additional, 
very red RC. This redward extension is an artefact of the 
observational dithering pattern in which some (less than 
1%) of the HB stars have poor i?-band photometry, hence 
this detection of very metal-rich, ancient stars is also an 
artefact, and they are excluded from our star formation 
history. 

Based on sp e ctrosc opy of 562 RGB stars in Fornax, 
IBattaglia et al.l (|2006f ) found two distinct populations: a 
metal-rich ([Fe/H] ^ —0.9) component, and a metal- 
poor ([Fe/H] ~ —1-7) component with a large metallicity 
spread. Consistent with observations of other dSphs, they 
found the metal-rich stars to be more centrally concen- 
trated than the metal-poor population. To examine the 
metallicities of the oldest stars, we show histograms of the 
metallicity distribution of the stars in the oldest age bin in 
Fig. [ini In Region 1, the metallicity distribution of these 
ancient stars peaks at [Fe/H] ~ — 1, but also contains many 
of stars with lower metallicities. The peak at [Fe/H] ~ — 1 
is also apparent in Region 2, but becomes less strong in 
Region 3 and is practically absent in Region 4. This outer 
bin appears to harbour a broader peak at [Fe/H] ~ —1.5, 
and a third peak at [Fe/H] ^ — 2. It should be noted that 
this third peak is most likely an accumulation of more 
metal-poor stars outside our metallicity range. However, 
the histograms of the inner three regions are consistent 
with the presence of three peaks at metallicities of [Fe/H] 
~ -1.0, -1.5, and < -2.0 dex. 
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Fig. 10. — Histograms showing the metaUicity distribution of 
stars in the oldest age bin. The histograms for the radial bins are 
labeled 1 through 4 going outwards from the center of Fornax. La- 
beled as "Total" is the sum of the four radial distributions. 

To summarise, we have recovered the metal-rich compo- 
nent (centred at [Fe/H] ~ —1.0) of Fornax, however our 
results support the presence of three distinct peaks in the 
metallicity distributi on of the ancient star s rather than 
the two proposed bv iBattaglia et all ()2006[ ). To test the 
statistical significance of the three peaks, we created ten 
Monte Carlo realisations of the Fornax stellar population 
using our SFH and age-metallicity relation. A comparison 
of the resulting metallicity distributions indicated that, 
despite some slight differences, the three peaked function 
was present in every synthetic populatio n. We therefore 
argue that the metal-poor component of IBattaglia et al.l 
can be sub-divided into two separate populations. 
This suggests that Fornax experienced three main star for- 
mation events in the period >10 Gyr ago, a hypothesis to 
be tested with further spectroscopic data. 

3.3. The Inner Shell 

In a previous, wide-area survey, man eTall (|200l 
noticed a shell-like feature near the cen ter of Fornax in 
our fi eld 15 (see Fig. [1]). Subsequently, lOlszewski et al.l 
(|2006l ) confirmed the presence of the feature based on deep 
photometry obtained with Magellan, and determined an 
age of 1.4 Gyr and a metallicity of [Fe/H] ^ —0.7 for its 
stars. Because of the relatively small number of stars in 
this overdensity and the large contamination by the overall 
Fornax stars, ob t aining a SFH is difficult, as also shown by 
lOlszewski et al.l (j2006t) . Therefore, we opt for constrain- 
ing the properties of the stars in the feature using a sin- 
gle component (SC) fitting strategy, described in detail in 
Ide Jong et al.l (|2008f ). In short, simple stellar population 
models with a narrow age and metallicity range are fit to 
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Fig. 11. — Background-subtracted Hess diagram and stellar pop- 
ulation constraints for the inner shell, (a) Hess diagram of the 
inner shell feature in grayscale, with darker pixels corresponding to 
a higher number of stars. The overplotted contours show the stel- 
lar density of the best-fit single-component stellar population, (b) 
Contours showing the regions of the age-metallicity plane where the 
single-component fits lie within 1, 2, and Scr from the best fit; the 
best-fit values are indicated with an asterisk. 

All stars in a 2.'5 wide elliptical annulus containing the 
inner shell were extracted from the survey. This was then 
divided into two regions: the 17° arc containing the shell, 
while the remainder of the annulus was used as a control 
field. The control field-subtracted Hess diagram of the 
shell is shown in panel (a) of Fig. [Til where the MSTO 
and RC of the shell population clearly stand out. Panel 
(b) of the same figure shows the areas in the age-metallicity 
plane that lie within 1, 2, and 3cr of the best-fit values, in- 
dicated with an asterisk. We find the age of the stars in 
the shell to be 1.6 ± 0.4 Gyr and the metallicity [Fe/H] 
= -O.gtoJ dex. Our SFH results indicate a very low SFR 
for this age and the average age-metallicity relation pre- 
dicts a higher metallicity of [Fe/H] ~ —0.5 for stars of the 
general Fornax population of the same age. Thus, these 
results seem consistent with the inte rpretation of the shell 
resul ting from an accretion event (jColeman et al.l l2004l 
|2005() . rather than being part of the underlying stellar 
populations. This accretion hypothesis will be discussed 
further in a later publication. 

3.4. Luminosity History 

With the SFH in hand, it is possible to construct the 
CMD of Fornax as it would have looked at some time in the 
past. In this way, the total luminosity of the system can 
be traced as a function of time. Combining the SFH fits 
of all fields, the overall SFR as function of age and metal- 
licity was used to construct artificial CMDs for Fornax at 
various points in the last 10 Gyr, assuming in all cases a 
binary fraction of 0.5 and a Salpeter IMF. By extending 
the CMDs far enough down the LF, the flux of all stars can 
be used to calculate the total luminosity. Although several 
systematic effects hamper a very precise measurement of 
My, the value we obtain for the present day. My — —13.0, 
is ver y close to the literature value of My = —13.1 (|Mateol 
I1998D . The evolution of the V-band luminosity of Fornax 
during the past 10 Gyr is shown in Fig.[T21 This figure is 
for illustrative purposes only, and any uncertainties given 
would be estimates at best. However, we can see that while 
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Fornax has experienced a general trend towards increas- 
ing brightness as more gas is converted to stars, there are 
significant variations caused by bursts of star formation. 
We find that the total luminosity of Fornax has a range of 
at least 1 magnitude over a Hubble time, hence at some 
point it was less than 50% of its current brightness. 
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Fig. 12. — A schematic diagram of the luminosity evolution of 
Fornax over the past 10 Gyr. 

4. DISCUSSION 

The mechanisms of star formation and chemical enrich- 
ment in low-mass galaxies are not well understood. Qual- 
itatively, the process is simple to envision. Initially, gas 
collapses at the centre of a dark halo to form the first gen- 
eration of stars. The heaviest of these quickly evolve and 
blowout material, injecting chemicals and energy into the 
surrounding gas cloud. This stellar feedback causes the en- 
riched gas cloud to expand and, if the potential well is deep 
enough, then collapse back to the centre of the dark halo 
to create a new generation of chemically enriched stars. 
Thus, a cyclical process is established in which multiple 
generations of stars form with a progressive increase in 
chemical abundance. 

However, there are still man y unresolved q uestio ns re- 
garding dSph star formation. 'Hel mi et al.l (j^OOG") have 
shown that, in contrast to the Galactic Halo, the dSphs are 
conspicuously lacking stars with [Fe/H] < —3.0. This sug- 
gests that the gas sourcing the oldest stars in dSphs was 
pre-enriched. Given that all dSphs are thought to contain 
some fract i on of ancient stars (i.e. with ages > 13 Gyr; 
iHeld et al.l l200d) the initial enrichment must have been 
an extremely rapid process, via some mechanism which is 
unclear. 

Furthermore, the SFHs for each dSph arc remarkably 
different. Some dSphs contain purely old stars (ages > 10 
Gyr) and are characterised by a simple SFH (e.g. Draco), 
whereas others are dominated by intermediate-age stars 
and have been able to maintain multiple epochs of star 
formation and chemical enrichment (e.g. Fornax). This is 



the leading question concerning star formation in dwarf 
galaxies: why do Draco and Forn ax reside at the ce ntres 
of dark halos with similar masses ([Walker et al.|[2007l ). yet 
they differ in brig htness by a factor of more than ten? 
iMaver et al.l (|2006D discussed this point, and described two 
possible scenarios: (i) Either Fornax and Draco initially 
contained the same amount of gas, hence local effects have 
allowed Fornax to produce stars with a greater efficiency, 
or, (ii) Fornax had access to a tenfold greater reservoir of 
gas, and Draco experienced a similar (but scaled down) 
initial star formation. We note that these are extreme 
sc enarios, and something in b etween is not precluded. 

iGrebel &: Gallaghed (|2004f) found that the reionization 
of the universe did not cause the expected reduction in 
dSph star formation, hence 'local effects' are thought to be 
the dominant factor producing the variety of dSph SFHs. 
Differences in the level of tidal distortion, mechanical feed- 
back and gas infall experienced by each system are cited, 
however t he precise nature of these local effects is uncer- 
tain (e. g. iDekel fc Silkl 119861: iMaver et al.l l2006D. In this 
regard, iFerrara fc Tolstovl (|2000^ have noted that a dSph's 
dark matter fraction will infiuence its SFH. External forces 
(e.g. tidal and ram pressure stripping) can remove blown- 
out gas from a dSph, however a massive dark halo will 
allow the satellite to retain its gas. 

4.1. Ancient Stars 

Fornax contains a large number of stars and has a com- 
plex star formation history, hence it is an ideal object to 
compare to simulations. Firstly, our best fits show that 
the ancient stars (age > 10 Gyr) in this system have a 
mean metallicity of [Fe/H] « —1.4 (Fig.[B]). This indicates 
that the first few Gyr contained an intense period of star 
formation and chemical enrichment. At some level, this en- 
richment appears to have occurred throughout the whole 
body of Fornax: Fig. [5] indicates that the oldest stars in 
every radial bin contain a number of [Fe/H] ^ —1.0 stars. 
However, we also note a metallicity gradient in this ancient 
population, such that the central stars display a mean iron 
abundance approximately 0.3 dex greater than those in the 
outer reg ions. This is consisten t with the spectroscopic 
results of iBattaglia et al.l (|2006[ ). The three peaks in the 
metallicity distribution function (Fig. [TU|) are possibly evi- 
dence for three main bursts for star formation in the early 
Universe. In summary, while our results indicate that the 
first few Gyrs saw a swift chemical enrichment process in 
Fornax, they also show that this enrichment was enhanced 
towards the centre of the system. 

In general, our results for the first few Gyr of Fornax 
are well reproduced by models. iMarcolini et al.l (|2006t ) 
constructed a 3D hydrodynamic model describing gas dy- 
namics and chemical enrichment in a dwarf galaxy includ- 
ing the contribution of supcrnovae. The metalicity dis- 
tribution function we find for the ancient Fornax stars is 
similar to that produced by the Marcolini et al. model. 
Furthermore, they find that stars located towards the cen- 
tre of a dSp h are the product of a more efficient chemical 
enrichment (IMarcolini et al. 20Q8 ), as is seen in Fornax 
(jBattagha et al.ll2006f) . iSalvadori et all (|2008il presented a 
semi-analytic cosmological model following star formation 
in a dSph galaxy in a Milky Way-type environment, pre- 
enriched to an abundance of [Fe/H] ~ —3. They showed 
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Gas expelled to the outer regions (or even the halo) of 
a satellite system is generally expected to be removed 
by r am pressure s tripping and tidal distortion. How- 
ever, iBlitz fc Robishaw ( 2000} ) presented evidence that Hi 
clouds exist in many Local Group dSphs, situated up to 
10 kpc (or, approximately three tidal radii in the case 
of Fornax) from the centre of the system. This is sup- 
port ed by the strong evidence for Hi assoc iated with Sculp- 
tor (iGarignan et all 1998 : Bouchard et al.. ,2003. ). Indeed, 
iBouchard et all (|2006^ discovered an Hi cloud located to 
the North of Fornax. Whether this is associated with the 
dSph is not certain (the radial velocity of Fornax is incon- 
venientl y close to that of the M ilky Way in this direction), 
however iBouchard et al.l (j2006[ ) state a min imum mass o f 
1.5 X 10^ Mq at the distance of Fornax. iMaveil (|2005D 
noted that an object following an orbit with a low ellip- 
ticity (such as Fornax; Piatek et al. 2007) will be better 
able to retain its gas. It therefore seems possible that the 
massive, extended dark halo of Fornax could have allowed 
gas in the outer regions to remain bound to the system for 
an extended period. 

An external influence on star formation in a satellite 
system are tidal in teractions with the larger host galaxy 
(jMaver et al.ll200(j ). Tidal forces experienced by a satellite 
syste m as it orbits its ho st can induce hurts of star forma- 
tion (Barton et al.' 200 0|) as the interstellar gas clouds are 
compressed (Mihos fc Hernauistlil996.). Usi ng images 
over a four year epoch, iPiatek et al.l (|2007l ) measured the 
proper motion of Fornax and derived an orbital period of 



that a dSph experiences intense star formation in the first 
few hundred Myr, with multiple cycles of star bursts fol- 
lowed by gas blowout and infall, with accompanying chem- 
ical enrichment. Their mctallicity distribution function 
and mean metallicity are roughly equivalent to those we 
measured for the old stars in Fornax. It is clear that simu- 
lations are able to accurately reproduce the first epoch of 
star formation in a dSph environment. 

4.2. Intermediate Age Stars 

ISalvadori et al.l (|2008( ) find a rapidly decreasing star for- 
mation rate with time, such that approximately 2.5 Gyr 
after virialization the rate has fallen well below lO""* Mq 
yr~^. This corresponds to an age of ~9 Gyr, or the sec- 
ond data point in our global SFH for Fornax (Fig. [H]), 
where we measure a star formation rate o f approximately 
3 X 10 "'^ Mq yr~^. Hence, although the ISalvadori et all 
()2008[ ) model provides an excellent reproduction of star 
formation in a dSph such as Sculptor, they note that a 
more complex SFH (such as that seen in Fornax) requires 
a different set of conditions. 

In this context, we find the star formation rate in Fornax 
to be an approximately constant value0 of 3 x lO^'^ Mq 
yr~^ in the period from 9 to 4 Gyr ago. This period also 
witnessed a slow, nionotonic increase in iron abundance. 
However, star formation in Fornax experienced a sudden 
increase approximately 3 to 4 Gyr ago, jumping threefold 
to ~10~^ Mq yr~^ with an accompying spike in chemical 
enrichme nt. These r e sults all agree with the Ca ii triplet 
results of lPont et al.l (|2004[ ). Our results also suggest that 
this epoch of star formation was relatively short-lived, last- 
ing 1 — 2 Gyr, and was confined to the central ~0.5rf (1500 
pc) of Fornax (Fig. [S]). 

There are a v ariety of explanations for this continued 
star formation. 'Salva dori et al.l (120081 ) state that a re- 
finement of the reionization criterion would allow massive 
dSphs with a lower initial gas-to-dark ratio in their mod- 
els, thereby leading to less efficient mechanical feedback 
and more regular star formation activity. This could ex- 
plain the steady star formation seen in the period from 
9 to 4 Gyr ago, however it cannot account for the sub- 
sequent burst. As an alternative scenario, Fornax may 
have experienced an injection of new gas to fuel this next 
generation of stars. We have already presented evidence 
for a merger in Fornax, proposing that a gas-rich dwarf 
galaxy merge d with this system to f uel st rong star forma- 
tion activity (jColeman et al.ll200l. I2005D . However, the 
timing is problematic: the original scenario requires that 
this merger occured approximately 2 Gyr ago, or 2 Gyr 
after the suddent burst of star formation found here. Fur- 
thermore, our chemical enrichment history shows that the 
new burst of stars was accompanied by a sudden increase 
in abundance, hence this would imply that, prior to star 
formation, the metal abundance of the gas was at least 
that of Fornax and therfore unlikely to be of foreign ori- 
gin. We would therefore argue that the 4 Gyr burst was 
fueled by gas originating in Fornax. 

This scenario requires gas blown away by star for- 
mation to have resided in the outer regions of Fornax 
for at least 5 Gyr before collapsing back to the dSph. 

^ One should always remember that each datapoint in our SFH represents an average measurement, and hence will contain hidden complexities. 
A useful time frame estimate for a single cycle of star formation and chemical enrichment is ~250 Myr. 



3.21^;^ Gyr with an eccentricity of e = O.lSi;;;^^. Tidal 
forces scale as hence this orbit implies that Fornax 

will experience a distortion force change of at least 50% 
as it moves from pericentre to apocentre. A pericentric 
passage approximately 4 Gyr ago is possible within the 
current solution (we are grateful to S. Piatek for sharing 
his orbital code), however the uncertainties are currently 
too large for a fair comparison between orbit and SFH. 



4.3. Recent Activity 

Finally, we examine the recent activity in Fornax. Fol- 
lowing the 4 Gyr burst, there was another star forming 
event 400 — 600 Myr ago, and a more recent event approx- 
imately 100 Myr ago (Fig. O. These 100 M y r old stars 
have been previousl y noted by IStetson et al.l (|l998t ) and 
ISaviane et al.l (|2000D . and are the youngest stars yet ob- 
served in a dSph. However, as a new result, we have also 
detected a small number of younger stars. Table [T] indi- 
cates that in the period 10 — 100 Myr ago, ^1500 Mq were 
formed in Fornax. This is minute compared to Fornax as a 
whole (see Fig. [7]), yet it is tentative evidence that Fornax 
may have been forming stars almost to the present day. 
We ignore the abundance measurements for this youngest 
population, given that they are based on only a few stars 
in each age bin and have significant errors. 
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Fig. 13.— Artificial CMD of the < 100 Myr stellar populations 
compared to that of the Fornax core region. The left panel shows a 
typical colour-magnitude distribution for the young stars in Fornax 
as predicted by our star formation history. The red points trace 
a 100 Myr population (the 7.95 - 8.10 log (yr) age bin in Table P 
while the blue points represent the <100 Myr stars. The right panel 
shows the CMD of the core region of Fornax. 

To demonstrate the difBculty in detecting this smah, 
age < 100 Myr population, we show an artificial CMD 
of the Fornax young stars in Fig. [T31 The left panel is 
a 'typical' predictior(f| for the young stars based on our 
SFH. Extinction, photometric errors and completeness for 
the survey have been included. Note that the distribu- 
tion of the blue points (age < 100 Myr) is almost identical 
to that of the red points (age = 100 Myr), thus these 
populations are degenerate and it is not possible to distin- 
guish between them based only on photometry. However, 
it shows that such a very young population could easily 
have evaded detection up to now. Given the small number 
of young stars predicted by our SFH, and the difficulty 
in separating them from slightly older stars, we classify 
this as tentative evidence for an ultra-young population 
in Fornax. Nonetheless, we expect very little foreground 
contamination blueward oi B — R = 0.4, hence the bright- 
est blue stars shown in the right panel would argue that 
Fornax did indeed experience star formation less than 100 
Myr ago. High resolution spectra of these young stars is 
required to accurately determine their age. Lithium is an 
important age diagnostic for late-type stars as it is eas- 
ily destroyed in stellar interiors, hence a follow-up survey 
could target the re sonance doublet of Li i at 6708 A (e.g. 
iMontes et al.l[200l[ l in these young stars to demarcate the 
age of the most recent star forming event in Fornax. 

5. SUMMARY AND CONCLUSIONS 

Based on two-filter photometry to a magnitude of i? 
23, we used a CMD-fitting technique to derive the star 



formation history for the Fornax dSph. All dSphs contain 
some number of ancient stars, however these systems are 
known to display a wide variety of SFHs. Fornax formed 
a significant number of its stars in the early Universe (age 
> 10 Gyr) and subsequently experienced a constant star 
formation rate. This behaviour can be reproduced by 
the simulations. However, in the period 3 — 4 Gyr ago, 
Fornax experienced a sudden burst of star formation, ap- 
proximately three times the rate for the previous 5 Gyr. 
The cause of this activity is unclear. The star formation 
rate has been decreasing ever since, with smaller events 
400 — 600 Myr and 100 Myr ago. We also find tentative ev- 
idence for a small number of stars (total mass ^ISOOM©) 
which have formed in the past 100 Myr. Fornax contains 
the most recent star formation activity of any Local Group 
dSph. 

Strong radial gradients in the SFH are also evi- 
dent. As noted by previous authors ([Stetson et al.l [19981 : 
ISaviane et al.|[2000l ). recent star formation has been con- 
centrated toward s the centre of Forna x. This trend is seen 
in other dSphs ([Harbeck et al.|[2001h . Furthermore, we 
have found that chemical enrichment was more efficient at 
the dSph's centre. Even the oldest stars display a metal- 
licity gradient, such that the inner stars have a mean iron 
abundance approximately 0.3 dex greater than the outer 
stars. Indeed, the first few Gyr in Fornax appear to have 
been a time of intense star formation and chemical enrich- 
ment: the age > 10 Gyr stars display a mean abundance of 
[Fe/H] — —1.4. The oldest stars in Fornax also show three 
peaks in the metallicity distribution, possibly evidence for 
three main bursts of star formation in the period >10 Gyr 
ago. We find the metallicity to have increased monoton- 
ically in the period 9 — 4 Gyr ago, and then experienced 
a sharp increase in conjunction with the intense burst of 
star formation described above. 

Thus, while the first few Gyr of star formation in For- 
nax can be reproduced by models, the cause of the burst 
4 Gyr ago is unclear. We have previously proposed that a 
gas-rich dwarf galaxy merged with Fornax approximately 
2 Gy r ago to produce sub-structure (Coleman et al. 2004 , 
120051 ). however we cannot reconcile the timing of this event 
with the observed peak in the SFH. Therefore, we suggest 
that gas enriched and blown out by earlier star forming 
events Hi clouds settled in the outer regions of Fornax, and 
then re-collapsed to fuel an intense period of star forma- 
tion at the centre of the dSph. This event may have been 
caused by tidal interactions with the Milky Way during a 
pericentric passage. 

The authors thank the anonymous referee for their com- 
ments which have improved the manuscript. The authors 
also thank S. Piatek and C. Pryor for their helpful advice 
regarding the orbit of Fornax in different potentials, and 
we acknowledge the use of their orbital algorithms. We 
are indebted to H.-W. Rix for his helpful discussions, and 
the MPIA observing team who obtained the MPG/ESO 
2.2m data described in this paper. 

Facilities: Max Plank: 2. 2m 



This prediction is susceptible to a variety of assumptions such as the IMF and binary fraction Here we have used the same recipe as described 
above; a Salpeter IMF with a binary fraction of 0.5. We constructed a large number of these artificial young populations and chose one which 
is most indicative of the class. 
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Table 1 

Fornax Star Formation History 



Age Bin 


SFR 


[Fe/HJ 


Mass 


(log-yr) 


(xlO--''M0 yr-i) 


(IO^Mq) 


7.05- 


7.20 


0.16 ±0.09 


-2.2 ±0.9 


0.73 ± 0.40 


7.20- 


7.35 


0.02 ±0.01 


-1.9 ±0.5 


0.16 ±0.08 


7.35- 


7.50 


0.09 ± 0.04 


-1.6 ±0.6 


0.81 ±0.29 


7.50- 


7.65 


0.0 




0.0 


7.65- 


7.80 


0.0 




0.0 


7.80- 


7.95 


0.0 




0.0 


7.95- 


8.10 


0.15 ±0.07 


-1.9 ± 1.1 


5.5 ±2.6 


8.10- 


8.25 


0.0 




0.0 


8.25- 


8.40 


0.04 ±0.03 


-1.2 ±0.9 


3.0 ± 1.9 


8.40- 


8.55 


0.11 ±0.04 


-0.6 ±0.7 


11 ±4 


8.55- 


8.70 


0.24 ±0.09 


-0.8 ±0.5 


35± 12 


8.70- 


8.85 


0.75 ±0.15 


-0.4 ±0.3 


154 ± 30 


8.85 - 


9.00 


0.14 ±0.06 


-0.2 ±0.6 


42± 15 


9.00- 


9.15 


0.12 ±0.06 


-0.9 ±0.9 


51 ±23 


9.15- 


9.30 


0.17 ±0.08 


-1.1±0.8 


98 ±46 


9.30- 


9.45 


1.50 ±0.63 


-0.5 ±0.6 


1240 ± 514 


9.45- 


9.60 


9.4 ±1.9 


-0.9 ±0.4 


10900 ± 2200 


9.60- 


9.75 


3.0 ± 1.5 


-1.1 ±0.3 


5000 ± 2470 


9.75- 


9.90 


3.6 ± 1.7 


-1.3 ±0.3 


8400 ± 3880 


9.90- 


10.05 


3.2 ±1.1 


-1.3 ±0.4 


10500 ± 3440 


10.05- 


- 10.20 


5.4 ±1.0 


-1.4 ±0.3 


25000 ± 4610 



